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1.0 SUMMARY 


The objective of the work discussed in this report was to develop a 
dynamic digital model for use in conducting stability studies of compression 
systems associated with turbofan (dual-spool) engines. This objective was 
accomplished through modification of the existing General £lectric Dynamic 
Digital Blade Row Compression Component Stability Model. The modi- 
fications included the addition of compression system dual path geometry, 
the implementation of the matching requirements between the two spools, 
and the considerations for radial redistribution of the flow ahead of the fan 
bypass splitter. 

The analytical equations on which the fan-core integrated model was 
based were first derived in a volume-averaged form for a three-dimensional 
flow. Then, they were reduced to a two-dimensional form for axisymmetric 
flow. The final equations satisfied the mass, axial momentum, radial 
momentum, and energy conservation equations for flow through a finite 
control volume. Use of the radial momentum equation in the model assured 
the radial equilibrium between the flows in the fan hub and tip regions. 

The relative total-pressure loss coefficient and deviation angle blade 
characteristics data for the fan huo, fan tip and core compressor blade rows 
wers calculated from TF34 fan and core performance test data. Verification 
of the calculated blade characteristics was conducted by reproducing 100%, 
*^5%, 90%, and 80% clean inlet corrected speed lines for the fan hub and the 
fan tip; and 104. 3%, 99. 8%, and 94.6% clean inlet corrected speed lines for 
the core compressor. Validity of the compression component models was 
demonstrated by comparing the component results with the test data for the 
point of instability and the location of the stage where instability is initiated. 

The fan and the core were then coupled together and the integrated 
system was run to instability by separately throttling the fan at 99. 1% and 
88. 3% corrected speeds, and the core along its 104. 3% and 99.8% corrected 
speed lines. The integrated model was then verified by comparing the 
model bypass ratio with that predicted by the steady-- te cycle deck, and 
by comparing the points predicted of instability with those obtained during 
test program. The fan speed/ccre speed relationship that was needed 
during this program was determined from steady-state cycle deck pre- 
dictions for the TF34 turbofan engine. 
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2.0 INTB ODUCTION 


One of the most significant parameters in the aerodynamic design and 
operation of gas turbine engines is surge margin. Aircraft engine designers 
are continually striving to provide adequate surge margin to ensure success- 
ful steady and transient engine operation. Operating exp*- ^ has shown 

that th^ surge margin of a compressor or engine can be ’ graa by a var- 
iety of factors. One of the significant factors is inlet d. .cortion both its 
steady- state and time-varying forms. 

Since it is not feasible to experimentally evaluate a compressor or 
engine for all possible distortion combinations, the need for a stability 
(surge) prediction technique has evolved whereby the loss of surge margir 
or performance can be determined analytically. Heretofore, the majority of 
techniques have been empirical correlations of past engine and compressor 
experience. However, in the recent past, GE has identified an analytical 
technique which is currently under development for analyzing the response of 
a compression component to either spatial or time-varying distortions. The 
GE technique is basically a one -dimensional explicit time-marching techni- 
que applied to a >lade row-to-blade row model, utilizing the conservation 
of mass, momentum, and energy relations. 

Under NASA Contract NASS- 18526, "J85 Digital Dynamic Distorted 
Compressor Model" (References 1 and 2), the General Electric Dynamic 
Digital Blade Row Compression Component Stability Model was expanded to 
provide a parallel-compressor analysis of circumferential total-pressure, 
total-temperature, and combined total-pressure /total-temperature dis- 
tortions, and induces circumferential redistribution of distorted flows in 
blade-free volumes. The model has demonstrated the ability to accurately 
represent the performance of the individual blade rows when quasi- steady- 
state throttling is performed along a speed line. Accurate prediction of 
instability at the experimentally determined s . rge line is shown by the 
growth of internally generated unsteady disturbance As part of the effort 
under NASA Contract NAS3- 19854, "Digital Surge Prediction Analysis" 
(Reference 3), it was shown that these disturbances represent aerodynamic 
instabilities ^ ither than numerical instabilities since a linearized stability 
analysis demonstrated that the real part of some of the eigenvalues of the 
Jacobian matrix of the system of equations changes sign from negative to 
positive at the point of aerodynamic instability. This change of sign is 
indicative of disturbances being amplified rather than attenuated. Under 
USAF-APL Contract F336 15-75-C -2029, "Digital Model Analysis of Planar 
Pressure Pulse Generator (P^G) and Two-Stage Fan Test Data" (Reference 
4), the response of a fan to planar waves was predicted and compared with 
existing test results. Analysis of the comparisons showed that the model 
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accurately reproduced the wave shape, amplitude, and phase-angle re- 
sponses up to a reduced frequency of 0. 2. Accurate levels of planar-wave 
distortion transfer coefficient were predicted, and the proper trends of 
sensitivity to surge with frequency were predicted. 

However, to date all efforts have been directed towards modelling only 
turbojet compression systems or the compression components of tvirbofan 
engines. Turbofan, multi-spool engines are characterized by one or two 
compression components on each spool. NASA-Lewis Research Center, in 
recognition of the need to develop stability models of more advanced tech- 
nology compression systems such as turbofan engines, initiated the effort 
discussed in this report. The objective of this program was to obtain a 
dynamic digital blade row model of the TF34 compression system. This 
compression system is composed of a single-stage fan without inlet guide 
vanes (IGV)on the low speed spool, and a fourteen- stage compressor on the 
high speed spool. The IGV and the first five stages of compressor are vari- 
able and are scheduled as functions of corrected speed. The primary de- 
velopment item in this effort was to develop a method for properly handling 
the matching of the two spools. This was accomplished by allowing for 
radial redistribution of the flow in the blade free volumes ahead of the lead- 
ing edge of the splitter and upstream of the fan. 

The existing General Electric Dynamic Digital Blade Row Compression 
Component Stability Model, now known as the Dynamic Analysis Program 
(DYNAP), was modified to permit simulation of the TF34 compression 
system. This was accomplished by modelling the flow path geome^- 
TF34 fan, compressor, and combustor and making an allowance tc with 
two spools. The fan speed and core speed relationship was define, J -"^ing the 
results from cycle deck studies. This relationship was used to assure that 
the proper speed match was maintained between spools as either the core or 
the fan bypass was throttled up a speed line. 

The single-stage fan hub, tip and mass-averaged maps were analyzed to 
produce steady-state blade row characteristics (relative total-pressure loss 
coefficients and deviation angles) which are representative of the hub and tip 
regions and the mass-averaged performance of the fan for corrected speeds 
of 80, 90, 95, and 100 percent. The core engine interstage data was stacked 
for the fourteen stages and analyzed to produce pitch-line blade row char- 
acteristics (loss coefficient and deviation angle) for the compressor at 
corrected speeds of approximately 95, 100, and 1C5 percent. 

The above described modifications to the computer model were verified 
by throttling the fan hub, the fan tip, and the core compressor components 
up their respective set of clean inlet speed lines to demonstrate speed line 
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and surge point replication. The stage initiating instability in the conapress- 
or was identified and compared with interstage data results. The compon- 
ents were then coupled together and the system was run to instability by 
separately throttling the fan and the core to their component surge lines. 

The bypass ratio n.igration compared favorably with the results of cycle 
deck studies. 
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3. 0 TF34 ENGINE COMPRESSION COMPONENT ANALYSIS 


Preparatory to developing a TF34 Compression System Model, tlie 
TF34 fan and compressor were individually modelled, the stability char- 
acteristics of each were determined, and the results were compared with 
test data. These efforts are documented in this section of the report. 

3. 1 FAN COMPONENT 


The TF34-GE-100 fan which is a single str.;e fan composed of a rotor 
and o.n outlet guide vane (OGV) was tested as a full scale component in the 
General Electric Large Fan Test Facility located in Lynn, Massachusetts. 
The splitter /gooseneck and the fan bypass flow duct were simulated. The 
bypass ratio was controlled through the use of both the main discharge valve 
and the bypass discharge valve. 

The fan was throttled at constant corrected speed to obtain the speed 
lines and surge point®. The results of these throttlings are shown in 
Figures 1, 2, and 3 eis the overall, tip, and hub maps, respectively for 80%, 
90%, 95%, and 100% corrected speeds. 

There are two items which should be noted with respect to these figures. 
The first is that the tip map is based upon the flow entering the bypass duct 
and the hub map is based upon the flow entering the core. The overall 
performance map is then obtained by mass weighting the hub and tip per- 
formance parameters. The second is that the so called surge line is defined 
by a stress boundary below approximately 102% corrected speed and by 
aerodynamic instability (surge or rotp.uing stall) above this speed. 

The fan was modelled using thirteen volumes (see Figure 4) as follows: 

5 inlet volumes 
1 rotor volume 
1 interaxial gap volume 
1 stator volume 
5 discharge volumes 

This representation and the methods used for the analyses conducted during 
this study are similar if not identical to those reported in Reference 1 and 
3. The length of the volumes was set so as not to exceed tVie length of the 
longest blade row, the rotor. The blade characteristics in terms of rotor 
relative total-pressure loss coefficient and deviation angle were determined 
for the four corrected speeds of interest using the fan performance infor- 
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Figure 3. TF34 Fan Hub Performance 










mation in conjunction with the fan geometry. In the absence of test data, 
the stator was assumed to be lossless. The model representations of the 
speed lines are shown in Figures 5, 6, and 7 as the dashed lines for the 
overall, tip, and hub performances, respectively. The agreement between 
the test and model speed lines is excellent especially in the low flow regions 
of the speed lines (near surge) with some small deviations noted in the high 
flow regions of the speed line. 

Stability predictions were made using the techniques of Reference 3 and 
the results are shown also on Figures 5, 6, and 7 by the open and solid 
circles. The open circles denote the region of aerodynamic stability (all 
roots of the characteristic equation of the Jacobian matrix have negative real 
parts) while the solid circles denote a region of aerodynamic instability (one 
or more roots of the characteristic equation of the Jacobian matrix have a 
positive real part). The stability predictions based upon the overall perfor- 
mance characteristics appear to match the test data quite well (Figure 5) near 
to where aerodynamic surge was observed. As expected in the region where 
the surge line is set by stress limitations, the erodynamic instability 
occurs at a lower flow. Examination of Figures 6 and 7, the tip and hub 
characteristics respectively, shows that the fan stability limits are set by 
the tip characteristics. That is, as the fan is throttled to surge, fan surge 
(stall) will originate at the tip of the fan. This finding agrees with test ex- 
perience which indicates that the fan is tip sensitive. 

3. 2 COMPRESSOR COMPONENT 

The TF34-GE-100 core was tested in the General Electric Multipurpose 
Test Facility located in Evendale. The test assembly included an exit flow 
duct to simulate the combustor flow path and an inlet duct faired into the 
gooseneck to simulate the core inlet flow path. The core compressor is a 
14 stage axial compressor with variable IGV and stators 1 to 5. The IGV 
and stators are ganged and controlled as a function of corrected speed vidth 
a Tt £5 bias. Although the TF34-GE-100 core compressor has stage 7 and 
10 bleeds, the test was performed with the bleed ports blanked off. Clean 
inlet compressor performance and stall mapping tests were conducted to 
obtain the 94.6%, 99. 8%, and 104. 3% corrected speed lines which are shown 
in Figure 8. 

The compressor model consists of forty-one volumes. There are 29 
volumes consisting of the IGV, rotors 1-14, and stators 1-14; and 12 free 
volumes consisting of two volumes ahead of the IGV and 10 volumes ahead 
of the turbine diaphram. The a.xial lengths of the free volumes are chosen 
not to exceed the length of the longest blade axial chord. Figure shows 
the schematic configuration of the model. 


10 




ajnsssjj 


11 


Figure 5, Model Representation of TF34 Fan Overall Performance 
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Figure 8. TF3^ Core Compressor Perlormance Map 









The blade characteristics, relative total pressure loss coefficient and 
deviation angle, \\ere determined from the stage characteristics, work and 
pressure coeff^'cients, provided by the Lynn Compressor Design personnel. 
Figure 10 illustrates ♦■he per.'‘'‘rmance of the compressor as calculated using 
the derived blade characteri;. ics in pol'ynomial form. Excellent compari- 
son with the demonstrated performance is observed at the 94.6%, 99. 8%, and 
104. 3% corrected speeds. 

Stability analyses were performed using the linearized stability techni- 
que of Reference 3. The results sho'wn in Figure 11 indicate very good agree- 
ment between the demonstrated and predicted stability points for 94. 6% and 
104. 3% corrected speeds. The predicted stability point at 99.6% corrected 
speed is 0.2 Kg /sec (0. 5 Ibm/sec) high in flow as compared to the experi- 
mental point. It is also interesting to note that the predicted stability point 
at 94. 6% corrected speed dees not occur at the peak pressure ratio, but 
rather on the positive slope portion of the sp' ed line. 

A time -dependent mod^l of the compressor was constructed using 
DYNAP for the purpose of conducting throttling -to- surge simulations. The 
compressor was represented in exactly the same manner as it was for the 
linearized analyses. The objective of these simulations was to determine 
what correlation, if any, existed between model observed surge line in- 
stabilities and test observed instabilities. The results are given in Table 1. 

Table 1. Instability Initiating Stage 


Test Results 


Corrected Speed 

Model Results 

Surge 

Rotating Stall 

94.6 

10 and 1 Simultaneously 

9-10 

2 

99. 8 

10 

11 

- 

104. 3 

1 and 2 followed by 11 

13-14 

4 


Examination of the throttling to surge simulations showed the anticipated 
instabilities (based upon the amplification function of Reference '' develop- 
ing in the back end of the compressor. Test data obtained from a TF34 core 
engine test using t-wo high-response wall-static pressure transducers mount- 
ed approximately 150° apart and aligned axially showed the surge pulse to 
develop in the back end of the compressor. The surprising result obtained 
from the simulations was the near simultaneous occurrence of instabilities 
in the front end of the compressor. Reexamination of the test results showed 
that rotating stall cells in the forward stages occurred at approximately the 
same time as the rear developing surge pulses. 
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Figure lO. Comparison of TFj4 Core Compressor Performance 
Demonstrated and Calculated, 
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These results are further verification of the capabilities of the General 
Electric Dynamic Analysis Program (DYNAP) especially its ability to pre- 
dict the regions in the compressor where instabilities will develop during 
throttling. However, since DYNAP is a quasi one-dimensional model it doas 
not give the nature of the instability, that is, surge or rotating stall. This 
determination is left to obtaining and examining test data or developing and 
obtaining; results from multi -dimeusicnal raodels. 
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4. 0 TF34 ENGINE COMPRESSION SYSTEM ANALYSIS 


In this section of the report, the analytical model for the computations 
of TF34 compression system flow dynamic behavior is discussed. The 
model is capable of simulating the interaction between the fan and the core 
compressor. Radial exchange between the fan hub and tip regions of the 
flow was considered to allow for the variation in the bypass flow ratio. This 
objective was accomplished by extending the one -dimensional time -dependent 
model, reported in Reference 1 to allow for radial flow redistribution. 

A brief description of the fan-compressor integrated compression 
system will be given in this section and will be followed by a generalized 
three-dimensional macroscropic balance flow analysis. The equations will 
then be reduced to describe a two-dimensional axisymmetric flow. 

4. 1 COMPRESSION SYSTEM MODEL DESCRIPTION 

The TF34 compression system has a single stage fan without an inlet 
guide vanes, and a fourteen- stage compressor. The IGV and the first five 
stages of the core comoressor are variable and are scheduled as functions 
of the corrected speed. In addition to the fan and compressor, the model 
geometry, as shown in Figure 12, includes the fan inlet duct, the bypass 
duct, the gooseneck, and the combustor. The purpose of including these 
duct volumes is to assure that realistic propagation paths for unsteady dis- 
turbances are simulated. 

The compression system was divided in sixty volumes, as shown in 
Figure 13. The fan inlet duct is represented by four volumes; the fan in- 
cluding the interaxial gap between the - or and the stator, and the fan exit 
ducting to the bypass splitter by eight %olumes: the bypass duct by four 
volumes; the gooseneck by five volumes: the compressor by twenty-nine 
volumes; the compressor discharge diffuser by four volumes; and the com- 
bustor by six volumes. All volumes ahead of the bypass splitter, including 
the fan rotor and stator, were divided into tip and hub volumes with a split 
stream-surface as shewn in Figure 13. The split stream- surface extends 
from the bypass splitter to the fan inlet and divides the entering flow into 
tip and hub flow rates in a ratio approximating the flow averaged bypass 
ratio. 

4.2 DEVELOPMENT OF QUASI-THREE DIMENSIONAL DYNAMIC MODEL 


The concepts used in the one-dimensional model of Reference 1 are ex- 
tended to include the compression system associated with a turbofan engine. 
The difficulties dealing with such systems do not only arise from having 
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dual compression units, bu. also from handling different bypass ratios. 

With different bypass ratios, the annulus stream surface separating core 
and fan bypass flows will migrate radially. Instead of dealing with variable 
surface boundary, a fixed split stream- surface is considered and radial 
flow is allowed to cross this artifical boundary. In the follo-'in-' paragraphs, 
the integral equations of motion will be written in general form and then will 
be simplified to a degree suitable for our present model. 

4. 2. 1 Governing Equations 

TVie equations of motion are derived for a general control volume, V, 
covering a segment of a v,oncentric ring. Figures 14a and 14b show the 
notation used to identify the location of the crntrol volu**'* ’nd the geometri- 
cal characteristics used in the equations. Variables with superscript, ”i", 
"j”, or "k” will denote the variables at axial, circumferential or radial 
faces of the control volume, respectively. 

Flow rates and velocities used in the equations of motion are depicted in 
Figures 15a to 15c. Figvire 15a illustrates the flows entering and leaving 
the control volvime surfaces, while flow velocities normal and tangential to 
the control surfaces are shown in Figures 15b and 15c, respectively. In the 
following derivations, the bar above variables indicates control-volume- 
averaged properties. 

Mass Conversation 


The integral equation representing the conservation of mass can be 
written as follows: 
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Where W , W^ , W , are the mass flows into the axially -facing , tan- 
gentially-facing and racially-facing surfaces, respectively. 
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Control Volume i, j, k 
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Let p represent the spatiaJly averaged density over the control voliume 
= ~ / P dV (3) 
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Application of tliis definition to the mass conservatio:^ equation results in 
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Momentum Conservation 


The integral equation representing the conservation of momentum in a 
given direction can be written as; 
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Application of Equation 5 to the axial direction and use of Figure 16a 
result in the following axial moment'um equation. 
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The axial momentvim of the fluid inside the control volume can be 
written as 
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b. Angular Momentum Balance 


Figure l6. A; al and Angular Momentum Balance 



where the bar indicates volume averaged quantities and is the control 
volume length in the axial direction. Therefore, the axial momentum 
equation can be written in the following form; 
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The principle of conservation of angular momentum is similar to the 
principle of the conservation of linear momentum. Application cf Equation ? 
to the blade-free control volvime shown in Figure 16b and manipulation of 
the different terms result in the following equation: 
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It should be noted that this equation is only valid for blade-free volumes 
since the tangential blade force is not included. 


Upon writing the momentum balance based upon Equation 5 in the radial 
direction for the blade -free elemental volume shown in Figure 17 and upon 
inregroting over the control volume, the following radial momentum equa- 
tion is obtained (see Appendix A): 
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where 0 is the angular extent of the control volume and bars indicate 
volume averaged quantities. For a complete annulus control volume, 0 is 
equal to 27T. The two last terms represent respectively the centrifugal 
and the pressure forces acting on the flow passing through the control 
volume. 

For a bladed control volume, this radial momentum equation is valid if 
the radial blade forces are negligible. 

Energy Conservation 

Energy conservation may be obtained in the same manner as the mass 
and moment\im balances. With thermodynamic relationships it can be 
written in terms of an entropy balance. For a control volume with zero 
heat transfer to or from the fluid, the entropy balance is given by: 
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Here the entropy is produced inside a control volume by any irreversible 
process, for example, flow undergoing total-pressure losses. Application 
of Equation 11 to a flow passing through the control volume shown in Figure 
15 permits the following equation to be obtained. 
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Here Sp is the total entropy generated within the control volume. 


( 12 ) 


The terms enclosed by braces, | | in Equations 4, 8, and 12 are 
additions to the one -dimensional forms of the continiiity, axial moment\im, 
and energy equations given in Reference 1 
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4.2.2 Force, Pressure, and Entropy Production Terms 


This set of equations (other than being applicable to three-dimensional 
flows without heat ‘■ransfer), describes the state of a flviid in motion. To 
solve these equations, it is necessary to supply the caloric and thermal 
equations of state and expressions foi' F 2 » ’^M» ^F* Equation 8 

represents the blade force acting upon the fluid and includes contribution 
from the energy-producing torque term and the drag force term. from 

Equation 8 represents the mean pressure over k and k+1 radial facing sur- 
faces of the control volxime and is calculated from a linear combination of 
the volume entranc and exit axial static pressures. Sp of Equation 12 is 
the term which represents the total rate of irreversible conversion of 
mechanical to internal energy and, in the case of the model, represents the 
entropy production due to blade row losses. Details of the derivations of 
these three terms may be found in Reference 1. 

4. 2. 3 Radial Eqxiilibrium - Two Dimensional Dynamic Model 

The primary objective of the present effort is to establish a dyaamic 
model capable of physically modelling the dynamic interactions between the 
fan and the compressor in a turbofan compression system. To allow for the 
varying bypass flow in, the splitter region, radial redistribution has to be 
considered. A two-dimensional version of the previously derived equations 
can be used to dynamics lly analyze a turbofan compressor system, assum- 
ing that there is no circumferential redistribution of the flovr. The two- 
dimensional governing equations with radial eqmlibrium are summarized 
below: 
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The flow in the fan duct ahead of the splitter is treated as a -dimen- 
sional flow and is obtained by solving the above set of equations. However, 
core and bypass duct flows are ar-sumed one-dimensional flows and are 
treated using the same set of equations after eliminating the terms enclosed 
by the braces | [ . 

Flow inside a fan rotor passage is in general a three -dim ensi-inal flow 
and is usually complicated by the presence of blade twist and the existence 
of secondary flows. Since the rotor radial blade force is not negligioi® and 
since it is difficult to consider the secondary flow effect, the model will lot 
include radial redistribution in the rotor voluine. Thus, the rotor is con- 
sidered not to alter the radial distribution of the flow passing through the 
rotor control volume. In general, this approach is suitable and commen- 
surate with the approach used in deriving the fan tip and hub blade charac- 
teristics where one -dimensional flow an ' -sis is assumed. 


In contrast, radial flow redistribution vi-ill be considered in the fan 
stator. However, it will be handled in a simplified manner. It will be 
as.,umed that only the centrifugal and the radial pressure forces acting on 
the flow will affect the radial distribution of the flow in the stator volume. 

4.2.4 Boundary Conditions 

Four boundary conditions are necessary and sufficient to obtain the 
solution of the set of first order equations given by Equations 13 through 16. 
Three of these conditions will difine three flow properties ac the iniax. 
station, while the last condition will define a property or a relation of pro- 
perties at the exit station. 

At the compression system inlet, the air flows axially and therefore at 
the inlet the radial velocity component is zero. Inlet total pressure and 
total temperature are specified and held constant. With these inlet 
conditions and together with the density value, obtained fi om the first vol- 
ume averaged density, all the other inlet flow properties can be calculated. 
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In a turbofan compression system, the flow entering the system 
branches into two flows, a core flow and a bypass flow. Hence, exit bound- 
ary conditions have to be specified for each flow. At each exit, one bound- 
ary condition is necessary and sufficient. Usually;^ it is either the exit 
static pressure (Pg) or the exit flow fxirction (W v Tg/Pg). The choice 
depends on whether the flow exits to a uniform static pressure such as 
through a fan bypass flow duct nozzle or to a choked flow such as a t\irbine 
nozzle. In the latter case specification of the evi; flow function is a good 
simulation. 

4. 3 CALCULATION TECHNIQUE 


The following subsections will outline the method of solution adopted to 
solve the governing equations, the computational steps performed to obtain 
the flow variables of the entire compression system and finally the process 
followed to obtain the initial, the time relauced steady- state, and the time 
dependent solutions. 

4. 3. 1 Method of Solution 

The explicit second order time marching technique developed in Re- 
ference 1, is used to solve Equations 13 through 16. The technique is based 
on a Taylor series ejipansion of the volume averaged flow properties at 
time ’’t" to predict the flow properties at time ”t + At". If X represents a 
volume averaged flow property, then the second order Taylor series can be 
written as. 


X (t +At) = X (t) 

0 t 


3t 



(17) 


This equation indicates that both first and second time derivatives of the 
flow properties are needed. First derivatives are evaluated directly from 
Equations 13 through 16; however, the second order derivatives are eval- 
uated by differentiating the same set of equations with respect to time. 

4. 3. Z Computational Procedure 

The present model of radial redistribution was developed as an exten 
sion of and is quite similar to the Dynamic Digital Blade-Row Compress n 
Ccmponent Stability Model discussed in Reference 1. Figure 18 illustrates 
the computational procedure in a block diagram format. 

Block lisa statement of the required dependent variable information, 
that is, volume -aver aged density, axial flow, radial flow, and entropy which 
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are available from either an initial guess, a steady- state solution, or the 
solution known from a previous time step of a time dependent analysis. 

From the volume averaged axial and radial flows, and with the proper 
interpolation scheme for station flow values of Wj. and time 

derivative of the volume averaged density can be obtained for all the vol- 
\imes, as indicated in Block II, 

In order to calculate the first time derivatives of the other variables, 
station values of all the variables have to be known. Block III of Figure 18 
illustrates the interpolation between the volume-averaged properties in 
order to obtain station-value properties. 

In the case of blade-free vol\imes, where no blade forces or entropy 
production takes place, it is only necessary to calculate station velocity com- 
ponents in order to evaluate the equations of change. As shown in the lower 
branch of Block IV, the assumption of constant absolute flow angle across 
the volume is made. Total pressure, total temperature, and other desir- 
able parameters are also calculated as this point. A special case of the 
blade free -volume calculations is the imposition of the boundary conditions. 
At the model inlet, zero radial flow, constant total pressure and total tem- 
perature as well as constant entropy are maintained. At the exit of a model 
of a turbofan compression system, either a specified flow function or a 
specified static pressure boundary condition is imposed. Change of the core 
flow frinction or the bypass duct exit static pressure with respect to time 
will in fact represent a simulation of the compressor or the fan throttling 
process, respectively. 

As indicated in the upper branch of Block IV, the presence of a bladed 
volume requires the net axial blade force, entropy production, and station 
velocitv components to be calculated. Calculation of the net axial force and 
entropy production terms require knowledge of the loss coefficient and 
deviation angle. This information is available as polynomial representa- 
tions which are functions of incidence angle. Stationary blade rows are 
assumed to be lossless with constant deviation angle. 

Once the flow conditions at the stations are completely described, var- 
ious quantities of interest can be calculated such as stage coefficients, 
diffusion factors, etc. With all the necessary quantities on the right-hand 
side of the macrobalances available, the first-time derivatives of the vol- 
ume averaged properties can be calctilated as indicated in Block V. 

After the procedure in Blocks III, IV, and V is carried out for the tip 
and hub volumes, as indicated by Block VI, the second time derivatives of 
the properties can be calculated (Block VII). Expressions for the second 
time derivatives of the volxime averaged quantities are obtained from 
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8. Dynamic Model Block Diagram. 
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differentiating the macrobalances with respect to time. Analytical express- 
ions for the time derivatives of the station properties can be evaluated by 
interpolating between volumes and through use of the macrobalances. The 
same calculations are carried out for both hub and tip volumes (Block VIH). 

The procedure for calculating station properties and evaluating first and 
second time derivatives of the volume -averaged properties can be carried 
out for any number and types of volumes (Block X) and is not dependent on 
the particular geometry being modelled. Once these calculations are carried 
out for all the volumes, the solution can be advanced to the next time step 
through use of the second-order Taylor's series approximation. Block IX. 

As specified in Block XI, the technique can be repeated for as many time 
steps as required by the event being simulated. 

4. 3. 3 Initial, Steady-State, and Time Dependent Dynamic Solutions 

A prime interest of this study is to determine the turbofan compression 
system operation under both steady- state and time varying conditions. The 
solution procedure adopted to achieve this goal and outlined in the previous 
section requires the knowledge of an initial distribution of the flow pro- 
perties along the compression system. This distribution can be obtained by 
solving the one -dimensional version of governing equations (Equations 13, 

14, and 16) with the left-hand- side time derivative set to zero, and with an 
assumed value of the bypass ratio. During these calculations radial flow 
across the splitting surface is not permitted. The use of the same govern- 
ing equations to develop initial conditions to the system is adopted to assure 
consistency between the initial conditions and the subsequent calculations. 
The initialization is achieved in a separate subroutine that was added to tha 
main computer program. 

Since the assumed splitting stream surface ahead of the bypass splitter 
does not represent in reality the actual stagnation stream surface which 
separates the core and the bypass flows, the initial steady-state solution 
obtained in the first phase of the calculations does not represent the actual 
steady- state solution of the flow in the compression system. A time relaxa- 
tion technique is followed to obtain the actual steady- state solution. The 
initial solution is relaxed in time while keeping the inlet and exit boundary 
conditions constant and allowing radial flow exchange between flows in the 
hub and tip regions of the fan duct. The equations used in this phase of the 
calculations are the two-dimensional time -dependent governing equations 
(Equations 13 through 16). The solution is assumed to have converged to the 
steady- state solution when the change in the axial flow rate at any volume 
from one time step to the next is less than 0. 0001 percent. 
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Time -dependent solutions of dynamic operations are obtained by per- 
turbing the system from the steady-state operation case. Again, the two- 
dimensional time -dependent equations are solved using the calculation pro- 
cedures of the previous section. Time- varying inlet or exit boundary con- 
ditions are included during this phase of calcxilations. The response of the 
system as a function of time can be achieved from the resvilts of these 
calculations. It should be kept in mind that the bypass-ratio is never speci- 
fied, but results from the radial redistribution that occxirs at any moment. 

The throttling process of a compression system is a typical dynamic 
operation where the flow through the system changes through variation of 
the exit area. If the rate of throttling along a speed line is extremely slow, 
the system will follow very closely the steady-state speed lines. 

4. 4 Clean Inlet Flow Steady-State and Throttling Simulations 

The two-dimensional time -dependent model established in Sections 4. 2 
and 4. 3 was used to determine the clean inlet steady- state flow and the clean 
inlet flow stability limit of the compression system of the TF34 turbofan 
engine. The steady-state solution, satisfying the radial equilibrium of the 
flow in the fan duct, was obtained using the time relaxation method. Throt- 
tling simulation was then conducted by throttling the fan bypass duct along 
two different fan speed lines and the core along two different core speed 
lines. Fan tip, fan hub, and compressor blade characteristics used for 
these computations are the same as those used for the component stability 
predictions of Section 3. 

The core /fan speed relationship used in this study was obtained from 
the steady-state performance cycle deck data. This relationship is pre- 
sented in Figure 19. The relation implies that the fan corrected speed is 
dependent only on the core corrected speed; however, the data obtained from 
the cycle deck indicated a slight variation in the fan speed with the fan bypass 
exit area. This variation is less than t 1% and thus, is neglected in the 
present model. 

The size of the time step has to be determined prior to the numerical 
solution of the time -dependent equations. A time step of 10”^ seconds was 
chosen for the present computations. This value satisfied the CFL (Courant, 
Friedrichs, and Lewy) stability condition and, therefore, it assures the 
numerical stability of the present explicit scheme. Also, it was found that 
this time step size is small enough to assure quasi- steady throttling during 
the system stability studies. 

From the model computations, the bypass ratio that satisfies the flow 
radial equilibrium in the fan duct ahead of the bypass splitter will be deter- 
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miiied. The validation of the model will be based on the comparison of the 
bypass ratio with the steady-state value obtained from the cycle deck data. 

The results that will be j. ”3ented in the following paragraphs are 
given for some of the fan duct volumes. Figure 20 illustrates the nomen- 
clature used for the volumes and the axial stations of the fm duct. Blade - 
free volumes are indicated by the capital letter "F", while the tan rotor 
and stator are indicated by the capital letters "R" and "S”, respectively. 

The cLxial stations are denoted by the numbers 1 to 13:1 to 5 for the duct 
upstream of the fan rotor, 6 to 9 for the interaxial gap between the rotor 
and the stator, and 10 to 13 for the duct behind the fan and ahead of the 
bypass splitter. 

4. 4. 1 Clean Inlet Steady-State Flow 

Time-relaxed steady-state calculations were conducted for two cases. 
The data for the first case are 104. 3% core corrected speed, 99. 1% fan 
corrected speed, 136. 08 Kg /sec (300 lb /sec) inlet mass flow rate, and an 
initially assumed value of 5. 302 for the flow bypass ratio. The data for the 
second case are 99. 8% core corrected speed, 88. 3% fan corrected speed, 

116. 35 Kg /sec (256. 5 Ib/sec) inlet mass flow rate, and an initially assumed 
flow bypass ratio of 5. 316. For these cases, the initial value of the flow 
bypass ratio was obtained from the cycle deck data. 

The steady- state solution was obtained by relaxing in time the initial 
flow distribution until radial equilibrium was achieved. Convergence was 
assumed when the variation of the volume averaged axial flow rate between 
two consecutive time steps was less than 0. 0001 percent for each volume. 
Figures 21 through 30 show the flow variable changes during the relaxation 
process. 

Figures 21 and 22 represent the relaxation of the fan hub and tip volume- 
averaged axial flow rates as functions of time. These results are for 99. 1% 
fan corrected speed and similar results for 88. 3% fan corrected speed are 
shown in Figures 23 and 24. It should be mentioned that the calculated 
results are plotted at intervals of 400 time steps. By examining these results, 
the convergence of solution occurs after approximately 4000 time steps. All 
the volumes of the hub region, initially started with the same volume -aver aged 
axial flow rates, but after convergence they ended up with different axial flow 
rates. The same is true for the tip volumes although the hub and tip axial 
flow rates are different. This is due to the fact that the assumed stream- 
surface, that separates the flow in the fan hub and tip regions of the fan, does 
not coincide with the actual stream- surface that separates the core and bypass 
flows. 
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Figure 20. Stations and Volumes Nomenclc ire 











Number of Time Steps (in liundreds) 






FiguT'! 2k. Helaxation of Volume-Averaged Axial Flow Rate 
as Function of Time [Fan Tip Region). 


The timo relaxation of the hub-to-tip radial flow for the different fan 
volumes is shown in Figures 25 and 26. Positive values indica*-e radial flow 
across the split stream-surace from the hub to the tip region and vice versa, 
negative values indicate ladial flow from the tip to the hub region. The 
j results show that the fow adjusts itself both behind and ahead of the fan during 

the relaxation ..rocess. The large radial flow rate value for volume F4 in 
comparison with that for volume FIO should not be misjudged. Actually, if 
the difference in the axial lengths between the two volunxes is considered, it 
is found that the radial flow flux (flow rate per unit area) at volume FIO, 
ahead of the bypass splitter, is larger than that at volume F4, ahead of the 
rotor. 

The time relaxatioiv o' static pressure at the stator inlet and at the 
splitter leading edge are shown in Figures 27 through 30. Tip and hub 
average pressure values are presented for both 99. 1% and 38.3% fan cor- 
rected speeds. It is clear that the steady-state pressure levels in all esses 
were achieved in less than 20 milliseconds. After that time very small 
adjustments in the pressures are neeced to allow other flow variables to 
adjust to thei. steady-state values. 

The difference between the hub and tip averaged pressures indic.tvi :he 
level of the radial force acting on the fluid due to the radial pressure .grad- 
ient. During the steady-state time relsucation, the radial force changes 
until radial equilibriunri is achieved thereafter it maintains a constant level. 

Figures 31 and 32 show the distribution of the volume averaged axial 
flow rate at 99. 1% fan corrected speed for the hub and tip regions, respec- 
tively. Figures 33 and 34 show similar plots for the 88.3% fan corrected 
speed Cijse. The distributions of the axial flow rate are given as a function of 
time and their shape depends largely on the geometry of the split stream- 
surface. 

The distribution of the radial flow rate at 99. 1% and 88. 3% fan corrected 
speeds are shown in Figures 35 and 36, re .pectively. The results indicate 
radial flow ahead as well as behind the fan rotor. 

Figure 37 shows the distribution of the hub and tip average static pres- 
sure at 99. 1% fan corrected spr.ed. Both the initial and the steady- state 
values are given. 

4. 4. 2 Fan Bypass Duct Flow Throttling 

Clean inlet bypass duct quasi-static throttling was conductea for 99. 1% 
and 88. 3% fan corrected speeds. Figures 38 and 39 show tne throttling lo^ii 
on the performance maps fer the fan tip and hub, respectively. The hub and 
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P'iyul't! 27. Average Static Pressure at Stator Inlet Plane 
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Kigure i2. Distribution of Volume Averaged Axial Flow Rate 
along the Fan Tii) Region. 
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Figure J3« Distribution of Volume-Averaged Axial Flow Rate 
along the Fan Hub Region. 











Fitjuro Jk. Diatribution of Vol uine-Averaged Axial Flow Rate 
along the Fan Tip Region, 









Figure 35» Distribution of Hub to J'ip Radial Flow Kate 
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Figure 36. Di stribut iom of Hub to Tip Radial Flow Rate 
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Figure 37* Average Static Pressure Distribution 













Fan Bypass Duct Throttling - Fan Tip Performance 
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Figurii 39. Fan Bypass Duct ThroltJiriQ - Fan Hub Performance 





tip corrected flow- rates are based on the duct total area and are calculated 
using the rotor hub and tip inlet flow rates, total temperatures and total 
pressures. Examination of these figures indicates that the fan tip corrected 
flow rate decreases during the throttling procedure while the fan hub corrected 
flow increases. 

Further investigation was conducted to examine the effects of the bypass 
throttling on the distribution of flow variables. Figures 40 through 42 show 
the results plotted at different time steps for the 99. 1% fan corrected speed 
case. Curves at time t = 0 give the steady-state distribution of flow variables. 
Curves at time t = 140 milliseconds give the distribution of the flow variables 
at the fan stability limit. 

Examination of Figure 40 indicates that the axial flow rate in the tip 
region decreases with time for all the fan stations. The results shown in 
the same figure also indicate that at time t = 0 an increase in the axial flow 
rates occurs between stations 4 and 5. On the other hand, a decrease is 
observed between the same stations as the fan is throttled toward instability. 
This indicates that a spatial deceleration of flow occurs in front of the rotor 
tip as the model nears the instability region. 

Figure 41 shows the distribution of the axial flow rate in the hub region. 
Since the flow rate at station 13 is the amount of flow entering the core, 
the results indicate that the core flow does not vary significantly during the 
bypass duct throttling. This is consistent with the constant pumping require- 
ments of the core compressor. However, it should be observed that the 
flow entering the rotor, in the hub region, increases during the bypass duct 
throttling in order to satisfy the radial equilibrium in the regions ahead and 
behind the fan rotor. 

The radial flow distribution during the bypass throttling is shown in 
Figure 42. It is noted that there is no significant adjustment in the radial 
flow at the interstage gap, between station 6 and Furthermore*, a major 
part of the radial redistribution takes place ahead of the fan rotcr. 

The flow bypass ratio, defined by the ratio of fan duct flow to core flow, 
was calculated at several points during the throttling procedure, for both 
99. 1% and 88. 3% fan corrected speeds. Figures 43 and 44 show a compari- 
son between these results with those obtained from the cycle deck data. A 
very good agreement is indicated and as a matter of fact, the average differ- 
ence is ..sss than 1%. This demonstrates the validity of the present model 
^nd its capability of handling compression systems with bypass ducts. 

To investigate the model capability to predict the stability of the inte- 
grated compression system, a detailec. throttling study was conducted along 
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Figure ^1. Effect of Bypass Duct Throttling on the Distribution 
of the Axial Flow along Uie Far. Hub Region. 
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F'igurn 42, Effect of Bypase Duct Tiirottling on the Distribution 
of the Hub to Tip Radial Flow. 





Figure h'). Variation of Bypass Flow Ratio during the Bypass Duct Throttling 




Figure kh. Valuation of Dynass Flow Ratio during Bypass Duct Throttling 





the 99. 1% fan speed line. The time variation of the fan tip and hub correct- 
ed mass flow rates are shown in Figures 45 and 46, respectively. It is 
noted that up to time t = 140 milliseconds, the change in the flow is gradual 
and well behaved (except for a few converging oscillations at the initial 
stages of throttling). However, further throttling beyond the experimentally 
determined stability line showed tendency toward rapid change in flow rates 
and finally solution failure. 

In an effort to establish stability criteria similar to that of Reference 1, 
the fan tip amplification function (ratio of the rotor -volume-flow time 
derivative to the bypass duct exit-volxime-flow time derivative) was investi- 
gated for irregular behavior during the throttling simulation. Figure 47 
represents the time history response of the fan tip amplification factor. As 
the bypass duct is throttled, the parameter exhibits a rapid increase in 
the region of flow instability. Since the level of the exit flow derivative is 
constrr’ned by the imposed exit boundary condition, the behavior of the flow 
derivat. ve ratio indicates that internal perturbations in the fan rotor tip 
volume are amplified significantly as the model nears the stability limit. 

Once the derivative ratio had reached a level greater than 2 in the rotor 
volume, the solution was unstable and termination of the throttling process 
wovild not prevent the solution from failing. Therefore, this stabilitv criteria 
was adopted for use in the model. 

Figure 48 presents the ratio of the exit axial flow to inlet axial flow of 
the fan rotor tip and hub volumes at a fan corrected speed of 99. 1%. As the 
fan is throttled toward instability, the data illustrates that the tip exhibits 
rapid mass storage and the hub exhibits rapid mass evacuation. This 
means that more flviid enters the rotor in the tip regions than leaves. Since 
flow instability is usually associated with flow blockage, this finding indi- 
cates that the flow instability is encountered in the tip region of the fan. 

Bypass duct throttling at a fan corrected speed of 88. 3% is illustrated in 
Figure 49. The well behaved corrected flow variation with time indicates 
that the flow for this speed line is stable. A detailed view of the throttling 
is shown in Figure 50. It is noted that the throttling has continued beyond 
the experimental stress limit and beyond the fan component stability limit 
obtained in Section 3. 1 of this report. To further check the absence of in- 
stability, throttling was discontinued and the time calculations were continued 
further at the two points indicated by 16000 and 20000 time steps. No erratic 
behavior was indicated and actually the solution settled at the point where the 
throttling was terminated. These results indicate that even if the fan com- 
ponent is unstable by itself, it can be stable in the integrated system. 
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Figure U5» Variation of Fan Tip Corrected Mass Flow Rate 
during Bypass Duct Throttling. 
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Figure ^*6. Variation of Fan Hub Corrected Mass Flow Rate 
during Bypass Duct Throttling. 
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figure ^O. Details ol Faii Tip 88.3% Speed Line 
(Bypass Duct Throttling) 



4. 4. 3 Cor « *. Compressor Flow Throttling 


Core throttling simulation of the TF34 integrated compression system 
wcxs conducted for 104. 3% and 99.8% core compressor corrected speeds. 
Figure 51 shows the throttling locii on the compressor performance map. 

The throttling line at the 99.8% corrected speed matches closely the com- 
pressor compo-'ent demonstrated performance, while that at the 104. 3% 
corrected speed deviates slightly from the demonstrated performance. In 
genera.’, these restdts show a good agreement with the compressor com- 
ponent performance. However, it should be mentioned that in the present 
r del, the specific heat ratio has been taken equal to 1.4 for all the cases. 
Tli.s conditio.i attributed to the slight deviation observed at the 104, 3%. 

The compressor stability limit at 104.3% and 99,8% corrected speeds, 
as predicted by the present integrated system model, are shown by the solid 
circles in Figure 51. The results indicate a good agreement with those 
obtained by compressor component model of Section 3. 2 and shown by the 
open circles in Figure 51. 

The V .10.1 ge in fan performance during core throttling of a turbofan 
engine depends mainly on the level of the bypass ratio. For a high bypass 
ratio, the core flow ru'ze is small in comparison to the total flow rate. In 
addition, throttling along a constant core speed line usually involves small 
variation in the core flow. Therefore, for the above two reasons it was 
expected that core throttling would not affect significantly th<* flew in the fan 
duct. Actually, in the present model of the TF34 high bypa ratio turbofan 
engine, it was observed that the adjustment of the flow occurs just ahead of 
the splitter. Moreover, it vras also noticed that the decrease in the total 
flow rate due to the core throttling is averaged out at the front face of the 
rotor resulting in insignificant changes in the performance of the fan hub and 
tip regions. In other w'ords, the core throttling ot the TF34 integrated com- 
pression system does not ch*nge significantly the performance point of the 
fan. 
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Core Throttling of the Integrated Compression System Model Core 



5. 0 CONCLUSIONS AND RECOMMENDATIONS 


The overall objective of the present effort was to conduct a computer 
study on the stability of the TF34 compression system. The study was per- 
formed in two steps. In the first step, the stability of the fan and the com- 
pressor components was examined using the linearized and the time-depen- 
dent one-dim'^nsional stability techniques. In the second step, the stability 
of the fan* :ore integrated compression system was investigated using a 
two-dimt isional stabiLiy technique developed to account for the radial 
equilibriu’^i of the flow ahead of the bypass splitter. 

D'oring the component study effort, the fan tip, the fan hub, and the 
core compressor were individually modelled and the clean-inlet flow map of 
each was accurately reproduced. Fan corrected speeds of 80%, 90%, 95%, 
and 100%, and core compressor corrected speeds of 94.6%, 99.8%, and 
104. 3% were simulated. The stability limits of each component were pre- 
dicted and the results matched the test data to an acceptable degree. Fur- 
thermore, these studies have indicated that the fan instability originated at 
the tip and that the core instability initiated in both the front and back stages. 
These findings agreed well with the experimental data and in the cas- :>f the 
core agreed with the fact that rotating stall occurred in the front stages at 
approximately the same time as surge initiated in the rear itages. 

During the integrated system study effort, the dynamic model for the 
stability studies was developed The model was based on solving the govern- 
ing equations in time to describe flow unsteadiness in the compression 
system. The governing equations were first derived in a three-dimensional 
form and were then reduced to a two-dimensional form that included the 
radial equilibrium equation. Based upon the time -relaxed steady- state 
results and the throttling resvilts, it was concluded that this time -dependent 
model correctly simulated the flow interaction between the core and the fan 
and showed the manner in which the Eow adjusts radially ahead of the bypass 
splitter in order to account for the bypass ratio variation. 

The time-relaxed steady-state solution was obtained for 99. 1% and 
88.3% fan corrected speeds with corresponding core speeds of 104. 3% and 
99. 8%, respectively. Even thovigh the fan was only divided into two regions - 
hub and tip - the f7ow distribution results indicated a reasonablv good repre- 
sentatit.i of the flow field. Further, the results indicated radial flow and 
radial pressure gradients behind as well as in front of the fan rotor. The 
radi?'' -iressure gradients showed strong dependency on the flow bypass ratio 
in El -ion to the dependency on the fan tip and the fan hub stage character- 
istics. 
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Fan and core throttling options were included in the integrated model 
for stability study purposes. With the use of both options, stability limits 
were predicted for fan corrected speed lines of 88. 3% and 99. 1%, and 
core corrected speed lines of 99. 8% and 104. 3%. Predicted stability 
limits in all cases, except at 88. 3% fan corrected sp^sed, compared well 
with those obtained by the component study and those observed d\iring com- 
ponent testing. At a corrected fan speed of 88. 3%, the fan duct throttling 
of the integrated system indicated that the flow through the fan was stable 
throTighout whole speed line data range. This result indicates that the fan 
at certain conditions can be stable in a fan-compressor integrated com- 
pression system when it would be unstable as a component under the same 
conditions. 

When the fan stability limit is reached, the flow instability was obser- 
ved to occur in the fan tip region. Near the stall limit, a spatial decelera- 
tion of flow was noticed in front of the rotor tip. This flow deceleration 
subsequently caused the fan instability. 

In addition to the stability stuoies, flow distributions in both the fan 
hub and the fan tip regions during the throttling process were examined. It 
was observed that fan bypass duct throttling affected both the fan hub and the 
fan tip flows. On the other hand, the core throttling showed very small 
effect on the flow distribution in the fan. This is attributed to the fact that 
the TF34 engine is a high bypass turbofan engine, and the changes in the core 
flow at constant core speed will not have significant effect on the fan total 
flow rate. 

It was also observed that the radial flow in the fan interstage gap was 
not affected daring the fan throttling. This was expected since the fan rotor/ 
stator gap is small and most of the flow adjustments take place either ahead 
the rotor or downstream of the stator and in front of the bypass splitter. 

A comparison was made between the model bypass ratio obtained 
during the model throttling process and the results obtained from the steady- 
state cycle deck. The comparison indicated close agreement for both fan 
corrected speed cases that were investigated. 

Based upon the aforementioned results, it was concluaed that the pre- 
sent integrated compression system model successfully simulated the flow 
interactions in the fan-core integrated compression system. Further, iv is 
concluded that the present model offers a unique method for the studying the 
stability of a turbofan compression system. 

The capability of the present model can i enhanced to increase further 
its value. The results of the analysis discussed in the previous paragraphs 
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indicate that it is possible to calculate some of the characteristics of the 
Qow radial redistribution ahead of the bypass splitter by dividing the fan 
flow into only hub and tip regions. However, if one is interested in more 
details of the flow field, it is recommended that the flow field be divided into 
multiple radial sections in the region of the fan. Also, the capability of this 
model can be extended to handle three-dimensional flows by combining it 
with the circumferential redistribution of Reference 2. These modifications 
and extensions allow the model to simulate more complex geometry engines, 
such as the variable cycle engine with inner and outer bypass duct, and to 
handle flows with general inlet distortions. 
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APPENDIX A 


RADIAL EQUILIBRIUM EQUATION 


Upon writing the radial momentum balance based on Equation 5 for 
the elemental control volume shown in Figure 17C, one obtains: 




(PC dV) 
r 


r r r r 


+ ( 5W^ C‘ - 

z r z r 


+ (W, C - W C ) 
*1 ^ *2 '2 


+ (w. c ^w. c ) 

®1 " l ~^ ®2 “ 2 ~ 


+ S 5F 


(Al) 


The radial momentum inside this elemental control voltime can be 
written as 




dV = p C 5V = P C (S A ) L 
r r r r 


dV = § W L 
r r 


(A2) 


where the bar indicates averaged quantity and L is the control volume 
cnaracteristic lengl . in the radial direction. 


The fourth term enclosed in brackets on the right hand side can be 
expressed in the following form; 


W„ C + w„ C = W, C 49 (A3) 

“l ^ »2 “2 ^ “ 


where and are the averaged value of the circumferential flow rate 
and velocity. It should be noted that this term represents the centrifugal 
force acting on the control volume. 
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If we neglect all radial forces other than pressure forces, then S 5F 
can be written using Figure 17a as: * 

2 5F = (L r\e - L 
r z z 




+ ;L L P,-^~+ LLP, 
r z 1 2 r z 2 




, k _k k+l_k+l. . An 
(r P - r P ) L \d 

z 


+ — (P, + P. ) L L 

2 12 r z 


If we write the average pressure as 


(p 4- p ) 

2 ' 1 2 ^ 


-i- (P'^ + 


and 


^ k+1 k 

L = r - r 
r 

then the pressure force can be w*'. ’en as 

I SF = (P^ - + r'^) L (A4) 

r 2 z 

Upon substitution of Equations A2, A3, A4 into Al, the following 
equation can be obtained '.-^r the elemental volume of Figure 17. 


3(5W ) 

T 

d t 


1 
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(5W^C^ . 

r r r r 


(5W^ c' - 

z r z r 


(W- c - w . c ) + (w c \d ) 

^ e u 


(P^ - P*''*'') t r‘‘) L \e 

Z z 


(A5) 
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Similar equations can be written for all the elemental control volumes 
of the control volume bounded by the plane** j and j + 1, and shown in Figure 
17a. Upon summation of all these equatiors, the following equation is 
obtained. 


Z W 

] 

Z t 




r r r r 


(W^C^ 

z r z r 


^ r d r ' 




+ (P^ - + r^) L © 

c z 


] 


(A6) 


where 0 is the angular extent of the control volume. 
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